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Abstract 

Recent developments in quantum chemistry, perturbative quantum 
field theory, statistical physics or stochastic differential equations require 
the introduction of new families of Feynman-type diagrams. These new 
families arise in various ways. In some generalizations of the classical 
diagrams, the notion of Feynman propagator is extended to generalized 
propagators connecting more than two vertices of the graphs. In some 
others (introduced in the present article), the diagrams, associated to 
noncommuting product of operators inherit from the noncommutativity 
of the products extra graphical properties. 

The purpose of the present article is to introduce a general way of 
dealing with such diagrams. We prove in particular a "universal" linked 
cluster theorem and introduce, in the process, a Feynman-type "diagram- 
matics" that allows to handle simultaneously nonlocal (Coulomb-type) 
interactions, the generalized diagrams arising from the study of interact- 
ing systems (such as the ones where the ground state is not the vacuum 
but e.g. a vacuum perturbed by a magnetic or electric field, by impu- 
rities...) or Wightman fields (that is, expectation values of products of 
interacting fields). Our diagrammatics seems to be the first attempt to 
encode in a unified algebraic framework such a wide variety of situations. 

In the process, we promote two ideas. First, Feynman-type diagram- 
matics belong mathematically to the theory of linear forms on combina- 
torial Hopf algebras. Second, linked cluster-type theorems rely ultimately 
on Mobius inversion on the partition lattice. The two theories should 
therefore be introduced and presented accordingly. 

Among others, our theorems encompass the usual versions of the the- 
orem (although very different in nature, from Goldstone diagrams in solid 
state physics to Feynman diagrams in QFT or probabilistic Wick theo- 
rems) . 
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1 Introduction 



The attempts to clarify the mathematical framework underlying quantum chem- 
istry, solid state physics, quantum field theories (QFT) and connected topics 
such as -for example, the grounds for the functional approaches, the principles 
underlying renormalization- as well as attempts to deepen our current under- 
standing of widely used techniques (effective Hamiltonians, adiabatic limits...) 
are often hindered by very basic questions regarding the underlying mathemat- 
ical methods. This is particularly obvious when it comes to developing new 
tools, see e.g. our studies of enhanced algorithms and formulas for the com- 
putation of eigenstates of Hamiltonians when the (unperturbed) ground state 
is degenerate [U [21 [3] , or of the combinatorics of one-particle- irreducibility for 
interacting systems 0]. 

The present article focusses on diagrammatics. We argue that Feynman-type 
diagrammatics belong mathematically to the theory of linear forms on combi- 
natorial Hopf algebras, which allows to generalize the theory to a much wider 
setting than the classical one. We cover the usual examples corresponding to 
vacuum expectations over commutative products of fields where the propagators 
are represented by edges -this includes for example the various Goldstone-type 
diagrams and the perturbative expansions parametrized by Feynman diagrams 
in quantum field theories. However we go beyond and cover the case of ex- 
pectations over general states, which requires the introduction of generalized 
diagrams jSJ H] ■ The same combinatorics happens to provide a pictorial descrip- 
tion of cumulants in probability. The Hopf algebraic approach also allows to 
study expectations of products of free fields (Wightman fields) and derive new 
Feynman diagram expansions in this setting. 

In a second part, we prove general linked cluster theorems, covering all these 
cases and making as explicit as possible the link between the elementary algebra 
underlying the theorems and the graphical content (that relies on connected- 
ness in the graph-theoretical sense). Notice that we avoid deliberately functional 
methods (see e.g.[6|): although very efficient to derive the classical QFT linked 
cluster theorem (using a generating function in terms of an external source), 
they lack the generality and simplicity of the combinatorial proof. Moreover, 
they cannot deal with noncommutative Feynman diagrams because functional 
derivatives commute. In the process, we promote another idea. Namely, linked 
cluster-type theorems rely ultimately on Mobius inversion on the partition lat- 
tice. 
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we aimed at developping a mathematical framework to deal with products of 
Wightman helds and their average values over general states, one of the prob- 
lems this article addresses. 

Notation 

We use the following convention: since various products will be defined along 
the article (such as * or 0), when we want to emphasize what product is used in 
an exponential, a logarithm or any other operation, we put the product symbol 
in exponent, so that log* emans that we compute the logarithm using the * 
product, ir 0n means that we compute the n-th power of x using the ©-product, 
and so on. 

2 Free combinatorial Hopf algebras 

The objects we will be interested in are combinatorial Hopf algebras in the sense 
of Joni and Rota [8], that is, bialgebras which coproduct is of combinatorial 
nature (obtained by "splitting" generating symbols according to combinatorial 
rules encoded by remarkable "section coefficients" -in high-energy physics, these 
coefficients correspond roughly to the symmetry factors of Feynman graphs). 
More specifically, we will be interested in families of Hopf algebras corresponding 
physically to bosonic or fermionic systems, to the usual algebraic structure of 
quantum fields (equipped with a commutative product such as the normal or 
time-ordered product) and to Wightman fields. 

Since our results are more general than what would be required by applica- 
tions to quantum systems, we state them in full generality and will show later 
how they specialize to particular physical systems or mathematical problems. 
Let X be a fixed ordered set, X — {x\, ...}. In most applications, X 

will be infinite and countable, so that the reader may think to X as the set of 
natural numbers. 

The notion of combinatorial Hopf algebra, goes back to [§]. The general 
notion is ill-defined in the litterature (there are many natural candidates, but 
at the moment no convincing general definition). We choose here a simple and 
relatively straightforward definition suited for our purposes that reflects some of 
the natural properties one expects from the notion when the underlying algebra 
is free commutative or free associative. 

Definition 1 We call free commutative (resp. free) combinatorial Hopf algebra 
any connected graded commutative (resp. associative) and cocommutative Hopf 
algebra H such that: 

• (Freeness) As an algebra, H is the algebra of polynomials (resp. of tensors) 
over a doubly indexed (finite or countable) set of formal, commuting (resp. 
noncommuting), variables <fii(xs) where S runs over finite subsets of X and 
i = I. ..rife, where k = \S\ and where the sequence of the n^, i 6 N is a 
fixed sequence of integers. 
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• (Equivariance) The structure maps are equivariant with respect to maps 
induced by substitutions in X. In other terms, any substitution (that 
is, any set automorphism) a induces a Hopf algebra automorphism of H 
which action on the generators is defined by o~(<f>i(xs)) '■= (f>i(x a /s\). 

The 4>i( x s) are called the generators of H, the (commutative or associative) 
monomials in the 4>i(xs) form a linear basis of H and are therefore called the 
basis elements. 

We do not look for the outmost generality, and many of our constructions 
and definitions can be extended in a fairly straightforward way to more gen- 
eral systems. For example, one might consider free partially commutative Hopf 
algebras which generators 4>i{xs) would satisfy partial commutation rules (e.g. 
all the 4>j{xi) would commute for a fixed i, but without <j)j{xi) and <fik{xi) com- 
muting for i ^ I, and so on). The so-constructed Hopf algebras can make sense 
in various applications and inherit all the properties of free commutative or free 
combinatorial Hopf algebras that are required for the forthcoming reasonings, 
we refer e.g. to [5] for details on the subject. 

Some further remarks are in order. Recall first that, by the Leray theorem 
(see e.g. Prop 4 in (TH]), any connected graded commutative Hopf algebra H 
is free commutative, so that the assumption that H is freely generated as a 
commutative algebra comes for free when one assumes that H is connected 
graded commutative. 

The key point that makes the Hopf algebra combinatorial (and, as we shall 
see, suited to Feynman-type graphical reasonings) is that we assume that a set of 
polynomial generators is fixed and behaves nicely with respect to substitutions 
in X. 

The particular case where rik = for k > 1 corresponds to the classical 
situation where the only propagators showing up in diagrams are the ones that 
describe the free propagation of a particle. 

Let us mention at last that a more pedantical definition of combinatorial 
Hopf algebras could be given in terms of vector species, following the approach 
to combinatorial Hopf algebraic structures in [TTJ [T2J [T3] . 

A very important consequence of the equivariance condition is the following. 

Lemma 2 For any subset S of X , let us write Hs for the subalgebra of H 
generated by the ^(xt), T C S. Then, Hg is a Hopf subalgebra of H . Moreover, 
any automorphism a of X induces an isomorphism of Hopf algebras from Hs 
to Ht> where T := o~(S). 

The second property is a direct consequence of the first one since a (by 
definition of the induced map) induces an isomorphism of algebras from Hs to 
Ht- The proof of the first assertion follows immediately from the equivariance 
condition. Indeed, notice that an element of H or of H ® H is invariant by 
any substitution that acts as the identity on a finite subset S of X if and only 
if it is a polynomial in the ^(xt), T C S (or, in H ® H, a sum of tensor 
products of such polynomials). Now, since, for T C S, ^(xt) is invariant by 
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any substitution that acts as the identity on S, the same property holds also 
true of A((j)i(xT)), from which the Lemma follows. 

3 Some remarkable Hopf algebras 

Our favorite examples in view of applications to quantum chemistry, QFT and 
solid state physics are simple ones. However the generality chosen (allowing 
for example tii ^ 0) is natural to handle nonlocal interaction terms in the La- 
grangians (think for example of the quantum chemistry approach with Coulomb 
interaction). Our general approach also paves the way to a unification of QFT 
techniques (Feynman-type diagrammatics), umbral calculus (duality and lin- 
ear forms on polynomial algebras) and combinatorics of (possibly ordered) set 
partitions. 

For notational simplicity, we treat only commutative algebras in the usual 
sense, that is we do not treat explicitly the fermionic case (Grassmann or exterior 
algebras/ Fermi statistics). However, as it is well known, there are no difficulties 
in switching from a bosonic to a fermionic framework -it just requires adding 
the right signs in the formulas, see [HI [15], but handling simultaneously the 
commutative and anticommutative case would have required introducing con- 
sistently signs in all our formulas. For notational simplicity, we decided to stick 
to the bosonic, commutative, case, and to let the interested reader adapt our 
results to the anticommutative setting. We use the langage of particle physics 
and quantum chemistry (so that the bosonic algebra coincides with the algebra 
of polynomials). 

Definition 3 (Bosonic algebra) We write for the algebra of polynomials 
over the set of (formal, commuting) variables 

<pi (xi ),..., 0i (x n ), <f>k(%i), — , 4>k{xn), ■■■ 

where Xi el. These algebras are naturally equipped with a coproduct A : B i — > 
B ® B that makes them Hopf algebras (that is, A is a map of algebras). The 
coproduct is defined on the generators 4>i{xi) by requiring them to be primitive, 
that is: 

A{(j)i{xi)) := faixi) <g> 1 + 1 <3 <f>i(xi) 
and extended multiplicatively to V (A(xy) = A(x)A(y)). 

The lower index k should be thought of as the number of quantum fields 
showing up (in a very broad sense), whereas the Xi should be thought of as 
points, momenta or more generally dummy integration variables. For later use, 
we allow k = oo (to encode the countable set of eigenstates of a given many-body 
Hamiltonian) and will write simply B for B^ when no confusion can arise. 

Definition 4 (Coulomb algebra) We write for the algebra of polynomi- 
als over the set of (formal, commuting) variables <pi(xi), I < k, <j}(x^j), where 
i =/= j G X . The coproduct is defined by requiring the generators to be primitive. 
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This algebra is used in non-relativistic many-body physics and quantum chem- 
istry. The Coulomb interaction desribes the force between a charge at point Xi 
and a charge at point Xj. These two points are linked by the interaction, and a 
specific variable <fi(xij) is used to describe this connection. 

Definition 5 (Tensor algebra) We write Tj? for the algebra of noncommu- 
tative polynomials over the set of variables 

4>i{xi), ipi(x n ), ipk(xi), —,ipk(x n ), ... 

Xi G X. The coproduct is defined by requiring the generators to be primitive. 

Various other free combinatorial Hopf algebras possibly noncommutative 
but fitting in the general framework of the present article are described (or 
follow from the results) in HH] ■ Let us quote the Malvenuto-Reuntenauer 
Hopf algebra and various Hopf algebras of tree-like structures. Although the 
decorated version of the Connes-Kreimer Hopf algebra of trees in [12] may have 
a use for QFT, a more promissing path in that direction is certainly provided by 
|17j , the belief of which we do share : "ultimately we think that combinatorics 
is the right approach to QFT and that a QFT should be thought of as the 
generating functional of a certain weighted species in the sense of [IB]"- We 
restrict however the examples in the present article to well-established domains 
of QFT and leave further extensions to future work. 

4 Graphication 

Let us start by recalling the construction underlying diagrammatic expansions 
and show how it can be extended naturally in a noncommutative/nonlocal set- 
ting. We call this process "graphication" , by analogy with the "arborification" 
process underlying tree expansions in analysis and dynamical systems [19| 120] . 
In all the article, as a tribute to the physical motivations the ground field is the 
field of complex numbers (although the results hold over an arbitrary field of 
characteristic zero). 

The reduced coproduct A : H i — > H ® H is defined by Vx 6 H,A(x) := 
A(x) — x (£) 1 — 1 ® x. The coproduct and reduced coproduct are coassociative 
in the sense that (A ® H) o A = (H ® A) o A (and similarly for A). The 
iterated coproduct and reduced coproduct maps from H to H® 1 are therefore 
well-defined and written AM, resp. A^. 

We assumed in the definition of combinatorial Hopf algebras that the co- 
product is cocommutative: with Sweedler's notation A(x) = xW ® x^ 2 \ this 
means that x 1 - 1 ' ® x^ = x^ ® x^K Many of our forthcoming results could be 
adapted to the noncocommutative setting. However, this hypothesis is particu- 
larly usefulf when it comes to graphical encodings. 

Definition 6 The graphication map Q is the map from H to Q)(H® n ) sym C 

n 
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0(.ff® n ) defined by: 



n 




n 



Here, the superscript sym in 0(#®™) s y m means that the elements in the 



image of Q are sums of symmetric tensor powers of elements of H. We will use 
the canonical isomorphism from covariants ®{H® n ) sym := H® n / ' S n (where 

n n 

S n stands for the symmetric group of order n) to invariants 



to represent invariants using the bar notation (that is, y\ ® ... <g> y n in Hf^ is 
written [yi|...|y„]). Notice that, by definition (since we deal with covariants), 
for any permutation a, [yi\...\y n ] = [y CT (i) |-|j/<r(n)]- 

For example, in the bosonic algebra, abbreviating cf)i to <fi (a notation we 
will use without further comments from now on): 



where I runs over all partitions I\ ]J ... JJ/fe, k — l...n of [n] :— {1, with 
inf{i G Ij} < inf{z G Ij+i}- Or (isolating some components in the expansion): 



The same formulas hold in the tensor algebra and in the Coulomb algebras. 
Notice however that, in the tensor algebra, products are noncommutative, so 
that the order of the products in the monomials does matter. For example, with 
self-explanatory notation: 

bPi {xi)\i>2 (xi )^{x 2 ) 2 ] ^ [^i{xi)\^{x2)^2{xx)^{x 2 )]. 

We call brackettings the terms such as [(j){xi) 2 (j){x 2 ) 2 \<j){xi) 2 <j){x 2 ) 2 ]. The 
length of a bracketting is the number of vertical bars | plus 1. The support of a 
bracketting T is the set of the Xi showing up in T. For example, 



n 




^i)^)..#„)) - 53 j] <t>{xi)\...\ j] <t>{ Xi )\ 



I ieli i<Elk 



g(<P( Xl ) 4 <f>(x 2 ) 4 ) = [^(arOV^) 4 ] +mxiMxi) 3 t(x 2 ) 4 } + ... 
+ I8&(x 1 ) 2 ( j ) (x 2 ) 2 \^x 1 ) 2 ^x 2 ) 2 ] + ... 



g(4>i(xi)(/) 2 (xi)(l) 3 (x 2 ) 2 ) = [(j>i{xi)(j) 2 (xi)(t>3{x 2 ) 2 ] + ... 
+ [Mxi)\Mxi)Mx2) 2 } + 2[^i(a;i)^3(a:2)|02(a:i)03(a:i)] + ••• 



sup([(f> 3 (x 1 )(j)4(x 2 ) 2 \(j) 1 (x 1 ) 2 (j) 2 (xs) 2 }) = {x 1 ,x 2 ,x$}. 
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For each Xi £ sup(T), we write di for the total degree of the (j)j(xi)s in T (so 
that for T as above, di = 3, d 2 = 2, dg = 2). For later use, we also introduce 
the product of two brackettings, which is simply the concatenation product: 

[wi|...|m„] • [i>i|..>„] := [u 1 \...\u n \vi\...\v n ]. 

These ideas and notation extend in a self-explanatory way to arbitrary com- 
binatorial Hopf algebras. The only change regards the definition of the support 
and degree: Xi is included in the support of a bracketting whenever there exists 
a pair (J, S) such that i £ S and (j)j(xs) shows up in V so that, for example, 

sup{[4> 3 ,{x 1 )(j)i{x2) 2 \(j)i (xi) 2 (/)2(x 8 ) 2 |05(a;i : 5 i 8ao)]) = {xi,x 2l x 5 ,xg,x w }. 

Similarly, the degree dg, of xg in this bracketting accounts for the ^5 (2^1,5,8,10) 
term and is therefore dg = 3. 

The coefficients in the right hand side of the equation defining the graphi- 
cation will be referred to as symmetry factors. They are closely related to the 
structure coefficients for the coproduct (in the basis provided by monomials in 
the chosen family of generators, e.g. the 4*i(xj) for the bosonic algebra -these 
coefficients were called section coefficients by Joni and Rota but encode also 
the symmetries showing up in the coproducts. Whereas, for the bosonic and 
Coulomb algebras, symmetry factors are defined without ambiguity, for other 
algebras (the tensor algebra or general free combinatorial Hopf algebras) a given 
bracketting may appear in the expansion of Q(M) for various monomials M in 
the generators (for example, in the tensor algebra, [ip(xi)\ip(x2)] appears in the 
expansion of Q(ip{x\)xj}{x2)) and of Q{il)(x2)'4>(xi))) . In general, we will there- 
fore write for the symmetry factor of a bracketting T in the expansion of 
Q{M) and will write simply sr in the particular case of the bosonic and Coulomb 
algebras (where a unique M exists giving rise to such a factor). 

5 Graphical representation 

Perturbative expansions in particle and solid-state physics are conveniently rep- 
resented by various families of diagrams. Feynman diagrams are the most pop- 
ular ones, but there are plenty of other families with construction rules often 
slightly different from the one underlying Feynman diagrams. Just to men- 
tion one interesting feature, Feynman diagrams are usually independent of the 
time-coordinate of the vertices (this is because of the definition of the so-called 
Feynman propagators) , whereas other families showing up in solid-state physics 
take into account causality systematically to construct their diagrams. These 
ideas are particularly well-explained in |21) . to which we refer, also for a com- 
prehensive treatment of the zoology of diagrammatic expansions. 

Defining a graphical representation associated to a free combinatorial Hopf 
algebra depends highly on the structure and particular features of the algebra. 
We define various such representations, by increasing order of complexity, fo- 
cussing only on the algebras (tensor, bosonic, Coulomb) we have chosen to inves- 
tigate in depth. The reader who needs to construct other taylor-made graphical 
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representations will be able to do so easily using our recipes (for notational sim- 
plicity and to avoid pointless pedantry, we omit to introduce the most general 
possible definitions since the process of designing them is straightforward once 
the leading principles are understood on some examples). 

Let us mention that, when two brackettings U = [ui \...\u n ] and V — [v± |...|f n ] 
have disjoint supports, their product corresponds to the disjoint union of the 
corresponding graphs (this is the usual product on Hopf algebras of Feynman 
graphs in QFT, see e.g. [H]). 

5.1 Commutative local case: bipartite graphs 

We focus in this section on the bosonic algebra. Each Feynman bracket (say 
T = [(pi(x 1 )(j) 1 (x2) 2 \(j)2(xi) 2 (t>i(x2) 2 (j)2(x3M3(xi)(t>3(x2)<t)3(x3)}, see fig. [J) can 
be represented uniquely by a bipartite (non planar) graph with unoriented col- 
ored edges (i.e. by a graph with 2-coloured vertices and colored edges) according 
to the following rule: 

1. For each Xi £ sup(T), recall that we write di for the total degree of the 
4>j(xi)s in r. Draw a x;-labelled black vertex with di outgoing colored 
edges (the colors being attributed according to the indices j of the c/>j(xi), 
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e.g. a 4-edge black vertex for x\ with one 1-colored edge (solid line), two 
2-colored edges (dotted line) and a 3-colored edge (dashed line). 

2. Running recursively from the left to the right of T, for each term inside 
brackets and bars (e.g. ()>i(xi)(fii(x2) 2 , then (^(zi) 2 ^^) 2 ^^), then 
03(^1)03(^2)03(^3)), select randomly according to the colors and powers 
showing up in the monomials outgoing edges of the corresponding vertices 
(e.g. select one 1-colored edge from the x\ vertex and two 1-colored edges 
from the X2 vertex). Connect these edges to a new white vertex (e.g. a 
new white vertex with 3 outgoing colored edges) . 

These edge-colored bipartite graphs will be called from now interaction 
graphs. See [5J 0] for applications. 

Some simplifications are possible in many cases of interest. For example, 
when the X{ s Eire dummy integration variables and can be exchanged freely 
in any computation (e.g. of scattering amplitudes), the labelling of the black 
vertices can be omitted. 

Another classical situation (encountered with scalar quantum field theories 
such as 3 or cj) 4 and in statistical physics) is the one where k = 1. In that 
case, there is just one possible color for the edges, so that the coloring of the 
edges can safely be omitted in the definition of the graphs. In QFT, when 
several fields coexist (e.g. in QED, where k = 2), instead of using colors, 
practitioners use often different representations for the edges that depend on the 
fields involved (typically, in QED, edges associated to electrons are plain lines, 
whereas edges associated to photon propagation are represented by a succession 
of small waves). Of course, the use of colors or the use of different shapes for 
the edges are strictly equivalent and a matter of taste and habits. 

Another classical simplification occurs when the only brackettings of interest 
for practical applications are the ones where the only monomials showing up in 
the brackettings are products of degree 2 of the form <l>i(xj)(j>i(xk)- In the 
corresponding graphs, the white vertices always have two outgoing edges with 
the same color, so that these vertices can be erased -what remains is a graph with 
only black vertices and colored edges (that correspond, physically, to particle 
types). These are the celebrated Feynman graphs that one encounters in QFT 
textbooks. 

From now on we will identify systematically brackettings and the correspond- 
ing interaction graphs. 

5.2 Commutative non local case: tripartite graphs 

By nonlocal case, we mean that some n,-, i > 1 may be different from zero. 
The canonical example we have in mind is the one of QED in the solid state 
picture, that is with instantaneous Coulombian interactions (in the particle 
physics picture the interactions are local and encoded by products of fields in 
the Lagrangian; this corresponds to the commutative local case). 

For simplicity (and in view of the most natural applications), we assume 
that the only brackettings of interest are those in which the monomials showing 
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3 X 3 



Figure 2: Graph of the bracketting [<M x i)^i( x 2) 2 |02(^i) 2 0i(£2) 2 02(£3)|03(:ei,2,3 



up are either monomials in the (f>i(xj), either a (f>i(xs), \S\ > 1 (in other terms, 
no nontrivial products involving a <j>i(xs) should appear in the bracketting). 
Each Feynman bracket (say for example 

r = [Mxi)Mx2) 2 \Mxi) 2 Mx2) 2 Mx3)\<t>3(xi,2,3)}) 

can be represented uniquely by a tripartite (non planar) graph with two kinds 
of unoriented edges according to the following rule (see Figure 2): 

1. We distinguish between the two kind of edges in the following way: the 
first type of edge is drawn as colored lines (the color is an index i with 
di 0); the second type has no color and is drawn as a sequence of waves 
(as the photon propagators in QED). We call these edges respectively plain 
and wavy edges. 

2. For each Xi £ sup(T), draw a x^-labellcd black vertex with di outgoing 
edges. Colors are attributed to the plain edges according to the indices j 
of the 4>j(xi), the other edges correspond to <f>j(xs)s and are wavy edges. 
For example, we draw a 4-edges black vertex for x\ with one 1-colored 
edge (represented by a solid line in Fig. 2), two 2-colored edges (dotted 
lines) and a wavy edge. 
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3. Running recursively from the left to the right of T, for each term in- 
side brackets and bars (e.g. (j)i(xi)<j)i(x2) 2 , then 4>2(xi) 2 (j)i(x2) 2 , then 
03 (2:1,2,3)), 

• If encountering a monomial involving cf>i(xj)s, proceed as in the com- 
mutative local case: select randomly according to the colors and 
powers showing up in the monomials outgoing edges of the corre- 
sponding vertices (e.g. select one 1-colored edge from the x\ vertex 
and two 1-colored edges from the x 2 vertex). Connect these edges to 
a new white vertex (e.g. a new white vertex with 3 outgoing colored 
edges). 

• If encountering a 4>j(xs), select randomly a wavy edge outgoing from 
the Xi, 1 6 S black vertex, connect all these edges to new j-labelled 
grey vertex. 

These tripartite graphs will be called from now nonlocal interaction graphs. 
The particular case of solid state physics QED enters in this framework. More- 
over, the general case we consider seems new and of interest, since it should 
allow to treat QED computations in a complex background (e.g. with non triv- 
ial vacua). These concrete applications (that originated this work, together with 
questions and remarks by R. Stora) are left for future work. 

Great simplifications occur in many simpler cases of interest. The first sit- 
uation encountered in pratice is the one where rii = n 2 = 1. This corresponds 
roughly to the case where one class of particles is present (say electrons, up to a 
change from the bosonic to the fermionic statistics) and where interactions are 
encoded by nonlocal terms (Coulomb interaction lines). In that case, there is 
just one possible color for the edges, so that the coloring of the plain edges can be 
omitted in the definition of the graphs. Besides, since n 2 = 1, the grey vertices 
can also be erased, so that one ends up with bipartite graphs with two types of 
edges (plain and wavy). If one assumes further that the only graphs of interest 
are those with two outgoing edges, these edges can be also erased. One ends up 
with the familiar Feynman diagrams with black vertices only and two types of 
edges corresponding to electron propagators and Coulombian interactions. 

5.3 Noncommutative local case 

We focus in this section on the tensor algebra. Each Feynman bracket (say T = 

[4>l(x 2 )4>l (Xl)lpl (x 2 )\lp2(Xl)lpl (x 2 )f 2 (x 3 )lpl {x 2 )i) 2 {xi)\lp 3 (x 3 )^ 3 (x 1 )lp 3 (x 2 )]) Can 

be represented uniquely by a bipartite (non planar) graph with unoriented col- 
ored and locally ordered edges (i.e. by a graph with 2-coloured vertices and 
colored, locally ordered edges) according to the following rule (the rule will 
make clear the meaning of "local order" that corresponds to an order on the 
edges reaching a white vertex) (see fig. [3]): 

1. Proceed first as in the commutative case: for each Xi £ sup(T), recall that 
we write di for the total degree of the <pj(xi)s in T. Draw a Xi-labelled 
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black vertex with di outgoing colored edges (the colors being attributed 
according to the indices j of the <f)j(xi), e.g. a 4-edges black vertex for 
x\ with one 1-colored edge (solid line), two 2-colored edges (dotted lines) 
and a 3-colored edge (dashed line). 

2. Running recursively from the left to the right of T, for each term inside 

brackets and bars (e.g. ipi(x 2 )ipi {xi)ipi (x 2 ), then ^2(^1)^1 (^2)^2(^3)^1 {x 2 )i> 2 (xi)—), 
select randomly according to the colors and powers showing up in the 
monomials outgoing edges of the corresponding vertices and order them 
according to the order of the appearance of the colors in the (noncommu- 
tative !) monomial (e.g. select a 2-colored edge from the x 2 vertex, label 
it a; a 1-colored edge from the x\ vertex, label it b; a 2-colored edge from 
the x 2 vertex, label it c). Connect these edges to a new white vertex (e.g. 
a new white vertex with 3 outgoing colored and ordered edges, where the 
order is defined by the labels). 

These edge-colored and locally ordered bipartite graphs will be called from 
now free interaction graphs. 

The usual simplifications are possible in many cases of interest, we do not 
detail them and simply mention that, when the only monomials showing up in 
brackcttings are of the form tpi(x m )'tl'i(xi) , the graphical rule amounts to con- 
sider "classical" Feynman graphs with labelled vertices and colored and directed 
propagators. 

The noncommutative nonlocal case could be treated similarly by mixing the 
conventions for the noncommutative local and commutative local cases. The 
exercice is left to the reader. 

6 Symmetry factors and connectedness 

This section is devoted to a technical but fundamental Lemma that connects the 
symmetry factors of graphs with the topological notion of connectedness. The 
lemma is the ground for the proof of linked cluster theorems and is particularly 
meaningful in the noncommutative case (Wightman fields). We assume in this 
section that the Hopf algebra is an arbitrary free or free commutative combi- 
natorial Hopf algebra which generators are primitive elements (this condition is 
satisfied by all the combinatorial Hopf algebras we have considered so far) . 

Let us introduce first a further notation. Let £ be a basis element (a commu- 
tative or noncommutative monomial in the generators) of a combinatorial Hopf 
algebra which support S decomposes into a disjoint union T]j V (the support 
is defined as for brackettings). We write xt and xy and call respectively T 
and ^-components of x the two basis elements (possibly equal to zero) defined 
by Xt ® xy :— (Pt <S> Py) ° A(x), where Py stands for the projection on Hy 
orthogonally to all the basis elements that do not belong to Hy . 

The reader can check that this definition amounts to the following: to get 
xt, replace, in the expansion of a; as a monomial, all the 4>i{xK), K C V by a 1 
and all the <j>i(xK),K nV r ^0,_ft'nT^0bya zero (and similarly for xy). 
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Lemma 7 Let ^ be a basis element and T = T± ■ r 2 a bracketting or, equiva- 
lently, an interaction graph (of any type) such that sup{T\) R sitp(r 2 ) = and 
s*/0 (recall that ■ stands for the product of brackettings). Topologically, this 
amounts to assume that T decomposes as a disjoint union of graphs. Then: 

• The basis elements ^sup^r^ an d ^su P (r 2 ) are non zero - 

• Moreover: s* = s *=" p(r i' s *=" p(r 2) 

In the commutative case, this last identity can be abbreviated to sr = Sr 1 Sr 2 - 

Indeed, since the 4>i{xx) are primitive, the coproduct of a product, say 
a\...a n , of <f>i(xi()s is the sum of the tensor products a^.-.a^ ®aj 1 ...a,j n _ k , where 
{ii, ik} is a (ordered) subset of [n] and {ji, j n -k} its (ordered) supplement. 
The hypothesis s* 7^ ensures that sup(ip) = sup(Ti) ]j sup(r 2 ), and that in 
tp, there is no factor ^(xk) with KCiV ^ $ and KC\T The first assertion 
follows. 

To prove the second identity, let us first notice that sf = s r sup(ri> sup(r2) . 
Indeed, let us use the same notation as previously and write 4 f = a\...a n . The 
coproduct and its iterations are constructed by extracting disjoint subsequences 
out of the ordered sequence of the a^s. On the other hand, the basis elements 
showing up in Ti and T 2 belong to disjoint sets -the relative ordering of the 
basis elements with support in sup(T\) and in swp(r 2 ) in the expansion of \I/ 
do therefore not matter, which proves the identity. 

We can now assume without restriction (because of the cocommutativity) 
that T = [S 1 ,...,S fc ,Ti,...,T,] with ri - [Si,...,S fe ] T 2 = [Ti,...,T,]. Besides, 
since A is an algebra map, we have: 

A[ fc+ 1(* 8up(ri )*.« P (r 2) ) = A^^ sup{ri) )A^^ sup{ r 2) ). 

The multiplicity jir '■= (k + Z)!s* of T in A^+^Cfy) is therefore obtained by 
summing the coefficients of the tensor products (X a n\ , X^+l)) m A^+'l (\I>), 
where a runs over S k +i/ 'Stab({X\, X k+ i)) and (Xi,..., X k+t ) = (Si, S k , T 1} Tj) 
(here, Stab((Xi, X k+ {)) stands for the stabilizer of (X\, X k+ i) in S k+ i). 
However, since the coproduct is cocommutative, these coefficients are all equal 
and one can restrict the computation to the tensor products (S^m, ...,Spr k -\, T^m, T a n\), 
where a and /3 run over permutations in S k / Stab{{S\, ...,S k )) and Si/ Stab((Ti, Ti)), 
and multiply the result by the number of fc-element subsets in [k + l] = {1, k+ 
I}. We get finally: 



l lT = \ k JATiMr 2 



that is, s* = s *°" p(r i ) s ^ s " p(r2) orj i n words (and slightly abusively), "the sym- 
metry factor of an interaction graph is the product of the symmetry factors of 
its connected components" . Notice that the property holds true also in nonlocal 
and/or noncommutative cases. 
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7 Amplitudes and Feynman rules 

A linear form on a combinatorial Hopf algebra is unital if p(l) = 1 and infinites- 
imal if p(l) = 0. 

Let us recall, for example, how linear forms on , X — {x±, ...,!„,...}, are 
usually constructed. Let T be an arbitrary finite sequence of integers. For any 
polynomial in k variables, say P(y x , y k ) = Pii,-,ikVi —Vk > we 

l<i 1 + ...+i k <n 

write P(T) for the polynomial Y[ P{4>i{ x t)i ■••> <f>k( x t)) G Bfe < One can think of 
P as the interacting part of a Lagrangian. Natural forms should then be thought 

4 

of as related physical amplitudes. For example, for k = 1, a typical P is A^j, 
corresponding to the 4> 4 theory (see [23] for details). The Green functions of 
this theory are computed via the formula 



where p denotes the vacuum expectation value of the time-ordered product of 
fields (see section fT0.il) . 

Let us treat now a radically different example to show the ubiquity of the 
approach. Let here the role of P be taken by Pj(i), the interacting Hamil- 
tonian of time-dependent perturbation theory (see e.g. our fTJ [2]): Hj(t) :— 
e iH ty e -iH t e -e\t\^ e ^ ^ e ^ j-^ e ig envec t rs of H with eigenval- 
ues Ai < A 2 < ... < A„ < .... We assume for simplicity that the ground 
state is non degenerate (Ai ^ A2), although the following reasoning holds in 
full generality due to fTJ [2] ■ The computation of the ground state of the per- 
turbed Hamiltonian Ho + V relies on the computation of the quantities such 
as: Y =< ei\Hi(ti)...H(t p )\ei >, where t > t\ > ... > t v > —00. We set 
k = 00 and (j) 2p (tj) := e " lX ^-i^ < e p \, (f> 2p +i{tj) := e +iX ^-^\\e p >. 
Then, P/(t) = J2i,j V i,j4>2j+i(t)(f>2i(t), where V itj :=< ej\V\&i >■ The unital 
form corresponding to the computation of Y is given simply by the form on 



where (0 iaJ (hj)\(j>i v+1 (*2j+i)) := ^-(hj^i+iihj+i) if *2j is even and i 2j+ i 
odd and zero else. The value of the form p on odd products is 0. Of course, this 
example is purely didactical and for such a computation the use of the formalism 
developped in the present article is largely pointless. It becomes useful when 
the situation gets more involved. Actually, the simple requirement of taking 
efficiently into account the divergences arising from the adiabatic expansion 
may involve advanced combinatorial techniques, see, besides the articles already 
quoted, our [3]. 

It is well-known that, in many situations, Green functions such as the ones of 
the (f> theory split into components parametrized by Feynman diagrams. This 




00 



p(cf)i (to)...<t>i 2k+1 (t2k+i)) ~ n Vi 



*2j + l,»2j 



(<Pi2j(t2j)\<Pi 2j+ i(t2j+l)), 



0<j<k 
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property also holds for more complex theories and is best explained through 
Hopf algebraic computations. Recall first that, since H is a Hopf algebra, the 
set H* of linear forms on H is equipped with an (associative; commutative if H 
is cocommutative) "convolution" product: 

Vp, pe H*,p* p(x) := p(x(r))n{x( 2 -)), 

where we used the Sweedler notation A(x) = xm ® £C(2)- Notice that, if p and 
p are infinitesimal forms, p * p{x) := p(xriy)p(xr2})> where we use the nota- 
tion A(x) = xj\\ ® X{2}- By standard graduation arguments, the convolution 
logarithm of a unital form p is a well-defined infinitesimal form t on P. We 
extend such a r to a linear form (still written r) on ~p® n (resp. Vf y " n ) by: 
t(xi ® ... ® af n ) := r(xi)...r(x n ) or r[xi|...|x„] := t(xi)...t(x„). 

Summing up, we get, for X an arbitrary monomial (basis element) in H, 
and since r is an infinitesimal form: 

Proposition 8 (Feynman diagrams/rules expansion) For an arbitrary uni- 
tal form on H , we have: 

P {x) = cxv* T {x) = Tog{x) 1 

or: 

P (X) = J2 s rT(T), 
r 

where T runs over all the brackettings (or interaction graphs) in the image of 
X by the graphication map. 

The map r acting on the Ts is called a Feynman rule. Applying Lemma [71 with 
H as in Sect. [6] we get immediately (with self-explanatory notations): 

Lemma 9 Assume that T = T% ■ T2, then: 

^r(r)= s ^ T (r 1 ) s ^r(r 2 ). 
8 The combinatorial linked cluster theorem 

The combinatorial linked cluster theorem expands a linear form on a combinato- 
rial Hopf algebra into "connected" parts closely related to the topological notion 
of connectedness. In this section, we show that this expansion is a very general 
phenomenon related to Mobius inversion in the partition lattice. Notations and 
conventions are in Sect. [6] 

We first recall some general facts on the partition lattice and Mobius inver- 
sion that are familiar in combinatorics but probably not well-known by practi- 
tioners of Feynman-type diagrammatics and linked cluster theorems (the par- 
ticular case of Mobius inversion for the partition lattice we are interested in here 
seems due to Schutzenberger, we refer to [24] for further details and references 
on the subject). 
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For an arbitrary set S, partitions t := {T\, T k } of S (that is: T\ JJ ■■■ \\T k — 
S, where {J stands for the disjoint union) are organized into a poset (partially 
ordered set, this poset is actually a lattice -two elements have a max and a min, 
this follows easily from the definition of the order). We write \t\ for the length 
of the partition (so that \t\ — k) and abbreviate the partitions of minimal and 
maximal length, respectively {S} and {{s}, s £ S} to 1 and 0. The subsets Ti 
are called the blocks, and the order is defined by refinement: for any partitions 
t and u, t < u if and only if each block of t is contained in a block of u. 

The functions f(x, y) on the partition lattice such that f(x, y) ^ only if 
x < y form the incidence algebra of the lattice. The (associative) product is 
defined by convolution: (/ * g)(x,y) := f{ x i z )9{ z -,y)- The identity of the 

x<z<y 

algebra is Kronecker delta function: 8(x,y) := 1 if x = y and := else. The 
zeta function y) of the lattice is defined to be equal to 1 if x < y and 
otherwise. The Mobius function /i(x,y) is defined to be the inverse of the zeta 
function for the convolution product. It can be computed explicitly: for x < t, 
where t = {Ti, T k }, we have: 

/i(i,t) = (-l)W + l*l(n 1 -l)L.(n fc -l)!, 

where 7ij is the number of blocks of x contained in T^. Using the identities 
\x * £(0, 1) = C * A*(0> 1) = <5(0> 1) = 0j we recover in particular the useful combi- 
natorial formulas: 

E E(- 1 ) |i|+|s| (^- 1 ) ! -(*fe-i)! = o 

0<fc<|S| t 

where t runs over the partitions of length k of S and ti stands for the number 
of elements in the i-th block T, of t and, with the same conventions on t, 

E Ec- 1 ) 1 ^ 1 ^!- 1 )^ w 

0<fc<|S| t 

The key application of these notions is to inclusion/exclusion computations 
in the partition lattice. Namely, for an arbitrary function h(x) on the lattice, let 
us set: h(y) = h(x). This formula defines uniquely h and, in the convolution 

x<y 

algebra: 

(h = h * C) ^ (h = h * m) 
so that h(y) = Y, h(x)fj,(x,y). 

x<y 

Let now x be a basis element of H with support written S, where H is as 
in Sect. [6] Let p be a unital form on H. Recall the notation xt denoting the 
"T-component" of x for an arbitrary T C S. For an arbitrary set partition 
t = {Ti, Tfc} of S, we extend p to a function on partitions of 5 and set: 

p x (t) := p(x Tl )...p(x Th ). 
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Recall also the decomposition p(x) — X)r s f T (r)- The Ls are represented by 
diagrams, among which some are connected. We set p C onn{x) := J2r c s f c T (Tc), 
where the r c run over connected diagrams. We can then apply the machinery 
of inclusion/exclusion to p x and define p x . The combinatorial linked cluster 
theorem relates p conn and p x : 

Theorem 10 (Combinatorial linked cluster theorem) We have, for an ar- 
bitrary basis element in a free or free commutative combinatorial Hopf algebra 
which generators are primitive elements: 

Pconn(x) =p x (l). 

We have indeed: 

;s*(i) = 5>*(i)MM) 

where t runs over the partitions of S and (with the usual notations) p(t, 1) = 
(-l)l*l +1 (|i| - 1)!. On the other hand, p x {t) = p(x Tl )-p(x Tk ) and p(x Tt ) = 
^Sp T *r(ri). Lemma [9] ensures that Sp . r = s r Tl ...s r Tk , and we get finally: 

Ti * l ■■■ fe l 

^(i)= E (-i) k+1 (k-iy. £ s x ri .,„. rk T(T 1 )-...-r(T k ), 

TilJ...LIT fc =s ri,...,r fc 

where r, is a graph showing up in the expansion of p(xr 1 )- 

Let now \1/ be an arbitrary graph with s|, ^ 0. The graph decomposes 
uniquely as a union of topologically disjoint graphs ]J ... JJ where n is 
the number of connected components of '5. We write Si for sup(^i) and = 
{Si, S n }. We have to show that the coefficient of t(^) in the right hand side 
of the previous equation is equal to Sp if n — 1 and to else. The first property 
is immediate, since if L is connected it appears only in the term associated to 
the trivial partition 1 of S and therefore with the coefficient (— l) 2 0!sp = sf. 

The second property is slightly less immediate, but follows from the princi- 
ples of Mobius inversion together with Lemma [SJ Notice first that L appears in 
the right hand side of the equation in association to all partitions of which 5$ 
is a refinement. The coefficient of L is therefore: 

[ J2 (-l) |t|+1 (l*l - l)!]s r 

s*<i<i 

which is zero as a consequence of the identity £ * p(Sy, 1) = 5(S^, 1) = in the 
partition lattice. 

9 The functional linked cluster theorem 

A linear form p on a free or free commutative combinatorial Hopf algebra H is 
called symmetric if it is invariant by a bijective relabelling of the variables xt -so 
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that p(<l)i(x2) s <i)3(x5) 2 <f>2(x9) 3 ) = p(0i(a;4) 8 03(^2) 2 '/>2(a;i) 3 ), and so on. When 
X is ordered, the form is called quasi- symmetric if it is invariant by a (strictly) 
increasing relabelling of the variables, so that e.g. p{4>i{x2) s ' ^{x^) 2 ^{x^) 3 ) = 
p(4>i(x4) 8 (j)3(x5) 2 (j)2(xs) 3 ), but is not necessarily equal to p(<?!>i(2;4) 8 </>3(:r2) 2 </>2(3;i) 3 )- 

Let us consider now a symmetric or quasi-symmetric unital form p on H, 
where H is as in Sect. [5] with an infinite ordered index set X. For notational 
simplicity, we will assume that X = N. That is, H is an algebra of polynomials 
(resp. of tensors) over a doubly indexed set of formal, commuting (resp. non- 
commuting), and primitive variables 4>i{. x s) where S runs over finite subsets of 
X and i = 1...7ife, where k — \S\ and where the sequence of the n^, k G N is a 
fixed sequence of integers. 

We first generalize the construction of the "interaction term" P in Sect. [Jj 
as follows. We let V — P(T)tcx be a family of elements of H such that P(T) 
is a polynomial (resp. a tensor) in the ipi(xs), S C T. 

Definition 11 We say that V is admissible if and only if 

1. For any order-preserving bijection <fi from T to R, 4>{P(T)) = P(R). 

2. For any T and any partition U\\V — T (where U and V inherit the 
natural order on T), we have: 

X)w(A)tr® {P b )v = P{U)®P(V), 

b 

where P(T) — ^2 b PbPb is the unique decomposition of P(T) as a linear 
combination of basis elements. 

We let the reader check that the map P constructed in Sect. [7] satisfies this 
requirement. 

The composition of p with P is a "scalar species": the value p(S) := poP(S) 
depends only on the number of elements in S (in the quasi-symmetric case, 
an increasing bijection induces the identity p(S) — p{T)) so that if we set 
/3( | S'l ) := ^syfj the scalar species p is entirely characterized by the formal power 
series 

K x ) ■^^P{n)x n . 

n 

This straightforward remark connects QFT and many-body theory with the 
various algebraic structures existing on the algebra of formal power series. Al- 
though apparently uselessly pedantical, it is actually useful to understand how 
these structures connect to the ones existing on scalar species. 

There exists a Hopf-like structure on linear combinations of finite sets (see 
e.g. [HI [16] for various developments of these ideas and the related notion of 
twisted algebras). The coproduct is defined by: 

8{S):= U®T 
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where ]J stands for the disjoint union, whereas the product is simply induced 
by the disjoint union of sets (the product of two overlapping sets is not defined). 
These maps induce a convolution product written (to distinguish it from the 
convolution product of forms on H) on scalar species: for a, /? two scalar species 
we get: 

a 0/3(5):= E a(U)0(T), 

U\JT=S 

or: a fi{x) = a(x)(3(x). 

Theorem 12 (General functional linked cluster theorem) We have, for 
any unital natural form p on H and admissible V: 



io g (p(x))=x;E-^ T ( r «) a;B 



where r := \og*(p), TJj runs over the connected Feynman diagrams with vertex 
set [n] and Sp c := X) -^fc s r c ' where P([n]) — X Abb is the decomposition of 

n b ™ b 

P([n]) as a linear combination of basis elements. 



Proof. We have: 



log(p(x)) = (log^)(x) = E (1 ° ge> ? (W) ^ 

n 

/iLI-1J^=H 



1 ^ (-l)_ fc + x 

k 4^ 



': IlII'- » l I I r, 

(-l) fc+1 ^ 



'■■ I: II II' " ' ' i 



where T/ i runs over the Feynman brackets in the expansion of p{h), ii '■= 

and TiTj denotes the concatenation of two brackets (so that e.g. [(j)(xi)\(j)(x5) 2 (j)(xs))[4 > (x2) 3 \4'( x i)) 

[^i)!^) 2 ^)!^) 3 !^!)]). 

Now, let r = Ti]J...]jrp be the (unique) decomposition of a Feynman 
diagram showing up in the expansion of /3([n]) into a product of connected (non 
empty) diagrams. According to Lemma [7] and since V is admissible, for any 
partition ^! JJ ... JJ ^4/ of [p], we have Sp 1 ... sp] — sji, where := JJ 1^ 

Al A ' jeAt 

and a.i is the number of vertices of T A t ■ 

The Theorem amounts then to the following properties: the coefficient of 
r(r) in \og(p(x)) is -^x n if p = 1 (that is if the graph is connected) and zero 
else. The first property (the connected case) is obvious from the expansion. 
Let us assume therefore that p > 1. The property follows once again from the 
general properties of the partition posets: the equation ([IJ concludes the proof. 
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10 Examples 

10.1 Quantum Field Theory 

In the quantum theory of the scalar field, the underlying Hopf algebra is the 
bosonic algebra B? , where k is the number of fields and the elements of X stand 
for dummy position or momentum variables. 

A typical example is the </> 4 (scalar) theory with k = 1 (we write simply (f> 
for 4>i ). The form p computes expectation values of time ordered products of 
free fields over the vacuum: 



p{<j){x x )...(t>{x k )) :=< Q\T( ( />(x 1 )..^(x k ))\0 > . 

Problems arise when some of the XiS coincide; these problems are the subject 
of the renormalization theory, we do not address them here. The physically 
interesting quantities are the interacting Green functions 

P U{x 1 )..4{x k )e l i't' i ^ dx ) 
G(xi,...,x n ) := — , '-. ( 2 ) 

pi e tj' 4> i {x)dx \ 

The key point is that p = exp* r, where r is zero if its argument has de- 
gree different from two and r((f)(x)(f)(y)^ :=< 0\T (cf>(x)(f>(y))\0 > is the Feynman 
propagator. The convolution logarithm can then be written as a sum of Feyn- 
man diagrams, where the lines represent Feynman propagators and the vertices 
represent spacetime points Xi. It can be checked that the standard Feynman 
rules of quantum field theory [6] are exactly recovered by the convolution ex- 
ponential |25j . The linked-cluster expansion provides a simple way to deal with 
the denominator of eq. © 126) . 



10.2 Cumulants 

Let X\, X n> ■■■ be a sequence of random variables. The underlying Hopf 
algebra for this example is once again the bosonic algebra B™ . The form p is 
defined by: 

p{cj3(x %1 )...(j){x lk )) := E[X h ...Xi k ]. 

This example enters the general commutative local case. 

When all the XiS are distinct, p{(j){xi 1 )...(j){xi h )) can be expanded as a sum 
parametrized by Feynman graphs which are disjoint unions of elementary graphs 
made of distinct black vertices, each joigned to a unique white vertex. The con- 
nected Feynman graphs appearing in the expansion correspond to the cumulants 
E c [Xi 1 ...Xi k \. The combinatorial linked cluster actually shows that this graph- 
ical expansion is equivalent to the classical identity: 

k 

E[X x ...X n ]=E c [X l ...X n ]+ ]T \{E c [X al ...X a] ], 

AiU...lM fc i=l 
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where A\ U ... U Ak runs over the proper partitions of [n] and Ai — {a\, a*.}. 

When the Xi are copies of a given random variable X and setting P([n]) := 
E[Xi...X n ] (which is admissible), we recover, using the functional linked cluster 
theorem: 

< e x > c -1 = log(E(e x )), 
with the convention < 1 > c = 1 and < X n > c := E c [Xi...X n }. 

10.3 Quantum field theory with initial correlation 

In solid state physics and quantum chemistry, the initial state is generally dif- 
ferent from the vacuum. The physically relevant form becomes (with the same 
notation as in the first example) 

p{<t>{xx)..4{xk)) ■=< $|Tfo(*i)...^(a; fc ))|$ >, 

where |$ > is a general state. It is also possible to consider a mixed state instead 
of the pure state |<I>). Except for very special cases (quasi-free states for bosonic 
fields and Slater determinants for fermionic fields) the convolution logarithm r 
of p is then more complicated than in the first example. In particular, t can 
be nonzero if its argument have degree different from two. In quantum optics, 
expansions in terms of r are known as cluster expansions and they lead to much 
better convergence properties [27] • For the fermionic fields, the convolution 
logarithms t are equivalent to the cumulants of the reduced density matrices, 
that are strongly advocated by Kutzelnigg and Mukherjee [28l [29) [30 l l31]. 

The diagrammatic expansions can then not be done anymore using Feynman 
diagrams constructed out of Feynman propagators: see e.g. our [1] and require 
the full apparatus of generalized Feynman diagrams for commutative local case. 

10.4 Non-Gaussian measures 

Perturbative expansions in statistical physics for measures of Gaussian type can 
be performed using the usual Feynman graphs of Sect. 110.11 This is because 
the Wick theorem applies. When dealing with arbitrary functional measures 
this is not the case any more: higher cumulants (i.e. higher truncated moments 
or truncated Schwinger functions) have to be taken into account. 

Feynman graphs and linked-cluster theorems have been developed by Djah 
and coll. [5] in this framework. They were extensively used in several problems 
of probability theory [32j [33l [34] . These Feynman diagrams are equivalent to 
those of a quantum field theory with initial correlation. 

10.5 Free probabilities 

Free probabilities deal with the noncommutative local case and study linear 
forms on the tensor algebra J-^. In general, the graphs required to study such 
forms are free interaction graphs. 
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In practice, the theory of free probabilities focus often on linear forms with 
particular properties. This allows for various simplifications and typical prop- 
erties as far as the corresponding cumulant expansions and their diagrammatic 
expansions are concerned. In particular, the Speicher's notion of free (or non- 
crossing) cumulant is obtained from the moment generating function by Mobius 
inversion with respect to the lattice of noncrossing partitions (and not with re- 
spect to the lattice of partitions), see e.g. [35] for further details and references 
on the subject. 

10.6 Truncated Wightman distributions 

In axiomatic quantum field theory, the form used in section 110.11 are replaced 

by 

W(x!,...,x n ) :=< ^\(j) H {xi)...(f)H{xk)\^ >, 

where the operator product is used instead of the time-ordered product, the 
fields are written in the Heisenberg picture and |$) is the ground state of the 
interacting system. Such functions are called Wightman distributions or correla- 
tion functions. The main difference with the standard case is that the Wightman 
distributions are not symmetric, the order of the arguments is fixed. Still, there 
is a perturbation theory of Wightman functions that leads to non-commutative 
Feynman diagrams 36 ]. Their combinatorics is the same as for the standard 
case [37j . Thus, the corresponding convolution logarithm t is again zero if its ar- 
gument is not of degree two, but now T(<t>{x),<t>(y)) = (0\<t>(x)<t>(y)\0) is different 
from T(cf)(y),<p(x)). 

However, it is also possible to work at the non-perturbative level and to 
define the form p{4>{xi) . . . 0(x„)) := W(x\, . . . , x n ). In that case, the convolu- 
tion logarithm r is generally not zero if its argument is not of degree two and 
r{(j){xi) . . . (f>(x n )) is now called a truncated Wightman distribution. 

The definition of truncated Wightman distributions was first given by Rudolf 
Haag in 1958 [55] . We follow (up to the order of the variables) the definition of 
Sandars' paper [315]: "For n > 1 we let V n denote the set of all partitions of the 
set {1, . . . , n} into pairwise disjoint subsets, which arc ordered from low to high. 
If r is an ordered set in the partition P € V n we write r S P and we denote 
the elements ofr byr(l) < ••• < ?"(|r|), where \r\ is the number of elements 
of r. The truncated n-point distributions u>^, n > 1 of a state ui are defined 
implicitly in terms of the n-point distributions 

U)(cp(xi),...,(p(x n )) = ^2 II ^M^U)' ■ ■ ■ > X r(\r\))- 

This is exactly the relation between p = exp* r and r for non-commuting vari- 
ables. 

From the physical point of view, the truncation procedure eliminates the 
contribution of the vacuum state as an intermediate state |40i p. 271]. Truncated 
distributions have many desirable properties. For instance, they decrease much 
faster than Wightman distributions at large space-like separation [41] . 
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10.7 Nonrelativistic systems with Coulomb interaction 



Let us neglect here the problem of dealing with the Fermi statistics (which 
amounts essentially to introducing the correct signs in the definition of the 
Hopf algebra structures, see e.g. [HI [T5] ). 

Let us consider n electrons in a quantum system of non-relativistic electrons 
with a Coulomb interaction in the external field generated by nuclei. This is 
the standard approach of quantum chemistry and solid-state physics. 

We assume that the non-interacting state can be described by a Slater deter- 
minant and that the particle-hole transformation was used to deal with occupied 
states. The form is defined as in section [10.11 and the Green functions are now 

p(ip(xi) . ..ip(x n )ipi(yi) . ..^{y n )e- u j 
87re |r-r'| 

The main difference with quantum field theory is that the interaction is not 
local. Still, in the linked-cluster expansion, we want to consider that the points 
r and r' in / are connected. In diagrammatic terms, r and r' are connected by 
a wavy line. We introduced the tripartite graphs to deal with this important 
case. 



G(xi, . . .,x n ,yi, 
where 



I = e 2 dtdrdr' 
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